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The gravitational polarizability properties of black holes are compared and contrasted with their 
electromagnetic polarizability properties. The "shape" or "height" multipolar Love numbers hi of 
a black hole are defined and computed. They are then compared to their electromagnetic analogs 
^EM rpj^g Love numbers hi give the height of the /-th multipolar "tidal bulge" raised on the horizon 
of a black hole by faraway masses. We also discuss the shape of the tidal bulge raised by a test mass 

- - - , m, in the limit where m gets very close to the horizon. 
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^ , I. INTRODUCTION 

, In Newtonian gravity, the quantitative theory of the "gravitational polarizability" of elastic, self-gravitating 
bodies was pioneered by Love who introduced two dimensionless measures of the response of an elastic body 
to an external tidal solicitation. To define them, let us first decompose the external tidal potential into multipolar 
^ components, say ?7ext = S = S Tir'^Pi (cos 9) , where the coefficient Ti measures the strength of the ^-th multipolar 
I ' component of the external tidal field. The first "Love number", hi, measures, essentially, the ratio between the l-th 
multipolar component of the distortion of the shape of the considered elastic body and T; , while the second "Love 
number" , fc; , measures the ratio between the Z-th order multipolc moment induced in the elastic body and Ti . 
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In the membrane approach to black holes (BH's) 0, [3, 0, H, @, BH's are treated as elastic objects, en- 
^ , dowed with usual physical properties. This raises the issue of defining and determining the BH analog of the Love 
' numbers hi and fc;. Some time ago, Suen made an attempt at defining and computing k2, i.e. the quadrupole 
] moment induced in a BH by an external quadrupolar tidal field. He surprisingly found that, with his definition of 
I the multipolc moments of a distorted BH, the tidally induced quadrupole moment was opposite to the externally 
applied quadrupolar tidal field (which would mean k2 < 0, in contrast with usual elastic bodies, which have fc; > 0). 
This unexpected result might well be due to the inappropriateness of the definition of induced multipolc moments 
adopted in Q. In two recent works, 0, it was found that if one defines the fc; Love number of a BH as being 
, the mathematical continuation, up to a star compactness formally equal to 1/2, of the (linearized) fc; Love number 
of a neutron star, the resulting ki vanishes for all Vs. However, there are subtleties inherent in any definition of the 
multipolc moments of BH's, so that there is currently no unambiguous determination of the ki Love number of BH's 
(see, e.g., [ll| and for discussions of the various ambiguities in the definition of the multipolc moments of BH's). 
^ I In the present paper, we shall not try to address these subtleties, and we shall, instead, focus on the computation of 
hi, i.e. on the quantitative measure of the tidal distortion of the shape of a BH. By contrast to ki, we shall sec that 
there is no ambiguity in the computation of the first Love number hi of a BH. 

Pioneering investigations of the tidal distortion of BH's were performed by Manasse [13 and Hartle [l^, [ij] (see also 
[istV Further investigations are due to D'Eath, [H, and, more recently, to Poisson and collaborators [Tsl. [Toj. 
Though these papers explicitly discussed the quadrupolar {I = 2) tidal distortion of BH's, they did not consider 
higher multipolar orders of tidal distortion, nor did they explicitly compute the value of the /12 Love number. In 
addition, there are sign errors in some formulas of p^ . which would affect the computation of /i2- 

The first purpose of this paper is to compute (in the linear approximation to tidal effects) the "shape distor- 
tion" Love numbers hi of a BH, for all values of / > 2. Let us mention that a recent investigation of the Love numbers 
of neutron stars [§| has found that in the formal limit where the compactness c^^ ~ GM/R of a neutron star tended 
to the compactness of a BH, c^^ = 1/2, the c-dependent hi Love numbers of neutron stars tended to the BH hi 
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values determined in this paper. The sequence of /ij's is a way of parameterizing the tidal effects due to a disturbing 
mass m located far away from the considered BH. We shall also be interested in studying the distortion of the shape 
of a BH when a disturbing test mass gets very near the surface of the BH. In this respect, we will find it useful to 
compare and contrast the gravitational polarizability of BH's to their electric polarizability properties. Let us recall 
that Hanni and Ruffini [20[ pioneered the study of the electric polarizability properties of BH's, and introduced 
the notion of the charge density, cr, induced on the BH horizon by an external charge o. [This concept of induced 
charge density was one of the origins of the membrane approach to BH physics 0, H, m, M, H, 0] •] Below, we shall 
use the concept of induced charge density to define electric analogs of the hi Love numbers, which we shall compare 
to the gravitational ones. Then, we shall also compare the evolution of the gravitational, or electric, induced effects 
as the external tidally-influencing (test) mass m, or charge q, approach the horizon. [Ref [2l| has considered the 
tidal distortion of the horizon of a BH by nearby (moving) masses, but has mainly used a "Rindler approximation" 
which replaces the BH horizon by a null hyperplane in a Minkowski background. In Section V, we shall discuss the 
connection of our results to those of [2l| .] 

The paper is organized as follows. In Section II, the notion of gravitational shape Love number of a BH is 
introduced. It is defined as the dimensionless ratio between the l-th multipolar component of the deformation of 
the horizon geometry of a BH and the l-th multipolar component of the external tidal potential generated by a 
static axisymmetric distribution of faraway masses. In Section HI, one introduces the notion of electromagnetic 
Love number by considering the influence of a distribution of faraway charges on the charge density induced on the 
surface of a BH. The case of a test charge approaching the BH horizon is examined in Section IV, while the case of 
an infinitesimal mass approaching the BH horizon is considered in Section V. The paper ends with a summary of our 
main results. 



II. MULTIPOLAR TIDAL DISTORTION OF THE "SHAPE" OF A BH UNDER THE INFLUENCE OF 

FARAWAY MASSES 



We wish to describe a physical situation where a BH of mass M is immersed in a generic, stationary, axisymmetric 
tidal field generated by faraway sources. For simplicity, we shall only consider here the static (nonrotating) case. This 
situation can be described by the Weyl class of static axisymmetric vacuum solutions of Einstein's equations, see, e.g., 
[2^ and references therein. The line element of a Weyl metric, in Weyl coordinates, reads 



-e^^dt^ + e^f''-^) (dp2 



Einstein's vacuum equations then imply the following equations for the fmictions tJj ~ ipip, z) and 7 — 7(p, z) 
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The linear differential equation (j2ap is identical to the 3-dimensional axisymmetric (^-independent) Laplace equation 
in cylindrical coordinates. Therefore, we can think of ip as being (minus) the Newtonian potential U = Gm/r 
generated by an ensemble of axisymmetric bodies: ip = —U . Without loss of generality, the potential U can be 
decomposed as U ~ Um + J/cxt, where Um refers to the BH of mass M that we consider, and Ucxt accounts for 
the external contribution(s). We recall that Um (for a Schwarzschild BH) corresponds, in Weyl coordinates, to the 
Newtonian potential generated by a "rod", along the z axis, of linear mass density p = 1/(2G) and of coordinate 
length Az = 2GM. The two remaining nonlinear field equations for the metric variable 7 can be solved by means 
of a line integral. The superposition of two or more axisymmetric bodies implies the presence of interaction terms in 
the function 7: 7 = i{Um) +7(C/oxt) +7int(C^M, C4xt), where 7int is bilinear (and nonlocal) in Um and C/oxt- We shall 
consider situations where the "central" BH of mass M is in "equilibrium" within C/oxt, i-C where the "elementary 
flatness" condition limp_»o 7 = is satisfied along the portions of the z axis that touch the BH horizon. [As is 
well-known, a nonvanishing 70 = limp^o 7 can be interpreted as the presence of a supporting strut or string] . This 
condition implies the constraint that C/cxt takes the same value at the north pole (z = +GM) than at the south pole 
{z = -GM) of the BH [H, [13, HI: 
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The notation u for the common value of [/ext at the north and south poles is introduced here for later convenience. 
Since the Newtonian potential U obeys a three-dimensional (j) independent Laplace equation, [/cxt can be decomposed 
in an axial-multipole expansion in spherical coordinates as 
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Here, {rw,dw,(t>) denote the spherical coordinates associated (as if one were in flat space) to the Weyl coordinates 
(p, z, 0), i.e. p = rw sm.{9w), z = rw cos{9w)- 

The coefficients Ti in Eq. ^ measure the strength of the l-th multipolar tidal field within which the considered 
BH is immersed. We are aware of the fact that the definition of the Tj's used here heavily relies on the quasilinear 
properties of the Weyl solutions. We think, however, that, in the linear approximation to tidal effects, i.e. in the 
formal limit where TJ —^ 0, for all Vs, the Tj's become unambiguously defined (and gauge-invariant, as they can be 
read off at r ^ GM). In general, ?7oxt (and therefore each J]) can be thought of as containing the superposition of 
the tidal fields due to a stationary ensemble of faraway (axisymmetrically distributed) masses. A particularly simple 
configuration is that where C/ext is generated by a single BH of mass m, located along the z axis at z = D ^ GM. 
[In this limit, the effect of any "string" supporting the perturbing mass m (and the BH) becomes negligible [HI]- As 
D ^ GM, there exists a wide region (with GM <C rw ^ D) where space is nearly flat and where we can expand 
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This shows that, in this case, the normalization of the tidal coefficients T; is such that 
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Let us also indicate (from [2^) the structure of the Weyl solution describing a "central" BH of mass M, 
distorted by several BH's of mass nii, described by rods of density 1/(2G) located in the intervals [bi — Grrii, bi 
along the z axis. The solution is described by 

i 

7 = 7a/ + Text = IM +'YlTn, + iMi, 
i i 

where (using henceforth G = 1 for simplicity) 
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and where each mi-dependent contribution reads (when suppressing the index i for readability) 
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with "fmirrij givcn by a similar expression (obtained by replacing M 
Here, 
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In Eq. (jSbp . C denotes an integration constant which must be chosen (after having satisfied the condition ^) so that 
linip^o 7 vanishes on the portions of the z axis which touch the "central" BH. When the condition ([3]) is satisfied, 



and the constant C in (jSbp is suitably chosen, it has been shown [23|, |2J, [25[ that the locus p = 0, 2 € [— Af, M], is 
a smooth (Killing) horizon. To see that its structure is that of a distorted BH, one would need to replace the Weyl 
coordinates (p, z, (j)) first by Schwarzschild-like coordinates (i?, 9s, (j)), such that 



p = \/R^ - 2MRsm9s, 
z ^ {R- M) cos 6*5, 



(12a) 
(12b) 



and then by (horizon-regular) Kruskal (or Eddington-Finkelstein) coordinates. For our present purpose, which is to 
read off the distorted geometry of the horizon R = 2M (p = 0), it will be enough to consider the geometry ^ in 
Schwarzschild-like coordinates, i.e. using the following consequences of Eqs. (|9a)) . ()9bp . (fTTj) and ([12]) 
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2M 
R 



R^dOl 



Here, as indicated in the particular case of ([8]), T/^oxt and 7oxt are defined such that -0 
where -^m and 7m are the undistorted "Schwarzschild" values ([9a|) and (j9bp . 
It is easily seen that Eqs. ([2]) imply that the following equality 



(13a) 
(13b) 

e-2'Aoxt^2gj^2^^^^2^ (13c) 

i'M + V'cxt, 7 = 7A/ + 7oxt, 



[7e: 
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const = 2u 



(14) 



holds on the horizon {p 0). Here, u is, as in Eq. ([3]), the common value of Ucxt = ^V'cxt at the north and south 
poles. Restricting the geometry (|13cp to f = to = const and R = Ro = const (and then taking the limit Rq 2M), 
and using Eq. (|14p . allows us to read off the geometry induced on the horizon: 
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dp^ + f3{p)d,p^ 



ds1i?=2A/ = (2Afe")^ 

where we have introduced the convenient variable 

/i = cos05=(^) , (-1<A^<+1), 

and where the fimction f3{p) describing the distorted horizon geometry is given by 

f3{p)^{l-p')p{p), 
(3{p) = e^^^^-'^'''^-'') . 



(15) 

(16) 

(17a) 
(17b) 



Here Ucxtip) ~ ^^'cxt|p=o is the value of Ucxt along the horizon, and u is, as in Eq. Q, the common value of Ucxt at 
the north and south poles. 

Note that the prefactor giving the overall length scale of the horizon geometry is not the mass parameter M entering 
the Weyl metric, but rather the "blueshifted" mass parameter 



Mu = Me" 



(18) 



While the mass parameter M is equal [25| to the Komar mass of the central BH, Eq. pSj) shows that the irreducible 
mass 0, [23 of the central BH is 
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where A denotes the area of the horizon. If we were to consider a physical process where our ensemble of gravi- 
tationally interacting masses would change their relative positions in an adiabatic manner, we should be careful to 
maintain all the irreducible masses (i.e. all the individual BH entropies) fixed during the process. This would mean 
that the Komar mass parameters M, rrii would need to change as the relative positions change, so as to compensate 
the variation of the blueshift factor e" in Eq. (jl8[) . However, as our purpose is to study the fractional distortion of 
the horizon geometry associated to external tidal fields, we will, in the following, factor out the overall scale factor 
{2MuY in Eq. (fT51) and focus on the shape of the /3(/i)-depcndent conformal geometry defined by the square bracket 
on the right-hand side of Eq. ([T5]). 



It can be easily proven that the Gauss curvature K — — ^ 

(fT5)) is given in terms of the second /z-derivative of the function /3(/i) by 
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{2MuYK = --f3 (m) 
As a check on Eq. (PO]) . one can verify the Gauss-Bonnet theorem: 
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of the distorted horizon geometry 
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When the horizon is regular, i.e. when Eq. ^ is satisfied so that (3^^ 
necessary for a horizon having the topology of a 2-sphere. 



/3(+l)+/?(-l) . 
/J'^^ = 1, Eq. ^ yields / / KdA = A-k as 



Let us now insert in Eq. (j20p our general parametrization ([4]) of an external tidal field by the sequence of 
tidal coefficients Tj. To do so, we must replace the Weyl-associated "spherical coordinates" rw, Ow entering Eq. ^ 
in terms of p, z, before taking the (singular) horizon limit p — > 0, z — s- Af/i which defines the horizon value of J7ext 
entering the definition of the function /3(/x), Eq. (|17ap . For instance, the / = 2 term in Eq. ^ would be rewritten as 



3 cos^ 9w 



1 1 



(22) 



whose horizon limit is simply = {MjiY. More generally, it is easy to see (using P;(l) = 1) that the horizon limit 
of r\^Pi{cos9w) is simply z' = (Mp)', so that 



C^Gxt(p) — t^oxtlho 



Y^TiM'y}. 



(23) 



Subtracting the value of C4xt| horizon at the north pole p = +1 finally yields the following explicit link between the 
Gauss curvature of the horizon and the sequence of external tidal coefficients T/: 



{2M^fim = -\^, [(1 - M')e^^' 



(24) 



At this stage, it is natural to decompose also the (scaled) Gauss curvature in multipolar components, say 

(2M„)2i^(p) = ^QPi(p). (25) 



Before deriving a formula giving the coefficients c; of the multipolar expansion of the Gauss curvature, let us note the 
two general, exact results 
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Eq. (|26ap is a restatement of the Gauss-Bonnet theorem checked above. As for the result (|26bp . it follows from the 
general structure (|20)) with a function /3(/x) satisfying (|17a[) and (jl7bp . Indeed, using Pi(/i) = p and the orthogonality 
of the Legendre polynomials, ci is given by 
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Using the expression (|17ap . where f3 takes the value 1 at both endpoints, on easily finds that ci necessarily vanishes. 

The orthogonality properties of the Legendre polynomials, J^-^ Pi{pL)Pi' (/i) = 2(5;;' /{2l + 1), then yield the following 
expression for the general q's as functions of the sequence of the T;'s: 

c; = ly^^P^i^)-^ [(1 -M^)e^S<' . (28) 

Note that q is a nonlinear function of all the T;''s, with /' being a priori unrelated to I. To better understand this 
function, we can however expand the exponential in Eq. (|28p in powers of the Tj's. In so doing, and in using the 
orthogonality of Pi{n) to all the powers /i' when I' < I, one finds that ci is the sum of a "principal" contribution 
oc T;, plus two types of extra contributions: (i) linear in T;/, which arise only with /' = / + 2, / + 4, / + 6 etc; and (ii) 
nonlinear in T;/, which can be quadratic cx T^, cubic cx T^, etc. Summarizing, we have the structure (for / > 2) 

ci = niTiM^ + ni,2Ti+2M'+^ + 7ii,4Ti+4M^+^ + ... + ^ n;;';"T;'r;" Af''+'" + ... , (29) 
where the n;'s, etc., denote some numerical coefficients. 

In the "linear" approximation to tidal effects (i.e. in the formal limit T;/ — > 0) we can neglect the nonlinear contribu- 
tions (X etc. In addition, if we focus, for simplicity, on the case where the (minimum) distance D to the external 
masses creating [/oxt becomes large, we find, in view of T; ~ 1/1)'+^, Eq. ([7]), that the contributions proportional to 
T;+2n are parametrically smaller than the "principal" contribution cx T; by a factor (Af/Z?)^" which tends to zero. In 
this double limit, we conclude that a kind of Hooke's law is valid in that the l-ih. multipolar component of the horizon 
curvature, q, is proportional to the /-th tidal strength coefficient T;, i.e. 



where 

(2/ + 1) '•+1 



ni 



ci ~ niTiM\ (30) 

p^d^,Pl{^i)^[{l-^,''){^,'-l)\ = 



^ ^ \i + m + 1) / d^iPl{^,)^J} ^{1 + 2)(z + 1)- ^ 



Finally, we conclude that, in the leading approximation, the Gauss curvature induced by a general external tidal 
potential reads (remembering the exact results of (|26ap and (j26bp ) 



{2M^YK{„) 1 + ^(/ + 2)[l + 1)^M'T;P;(m). 



(20! 



(32) 



Let us now convert this result in terms of the "shape" Love number hi. We define this number, in the present general 
relativistic setting, by paralleling the definition Love used in Newtonian gravity The basic idea is to write, for each 
multipole I, the specific potential energy gQ{5R)i associated to the l-ih "tidal bulge" of height {5R)i as /i; times the 
external tidal potential, evaluated (as it would be in Newtonian gravity) at the undisturbed radius of the considered 
object , i.e. 

goiSR)i^hiUrHRo). (33) 

or, equivalently, 

Ro J I goRo 

Here, go = GM/Rq denotes the (Newtonian) surface gravity of the considered object. We recall that the same 
(Newtonian-looking) formula is also used in BH theory to define the (renormalized) surface gravity of a Schwarzschild 
BH, if, as we shall do, one defines the "radius" of a BH as its (areal) radius Rq = 2GM. [In view of the remarks 
above, and of our focus on the linear approximation to tidal efi^ects, we shall not worry here about the possible 
distinction between the Komar mass M and the irreducible mass M„ = Me" ~ AI{1 + u + 0(u^)).]. As for [/ext(-Ro) 
we define it, as it would be done in Newtonian gravity, by taking the analytic continuation down to R = Rq of the l-th 
multipolar order asymptotic Newtonian potential. In the (intermediate) asymptotic domain Rq <ti R 'ti D the external 



7 



tidal potential ([4]) can be written in Schwarzschild-type coordinates {R,9s), Eqs. (fT2|) . as Uoxt — X]/ '^'-^'^'(^os 6'5). 
Formally continuing this Newtonian- like expression down to R = Rq, with fixed Schwarzschild-likc colatitude Og, 
leads to 

C/ext(i?o, - E t^r'(^o)Pi(A^) = TiRoPiiy). (35) 

where, in agreement with Eq. (jlSp . we have identified cos(?5 with the variable ^ used in the study of the horizon 
geometry above. 

Summarizing, the "shape" (or "height" ) Love numbers hi are defined by writing 

5R\ TiR\^ TiR^^^ 

In order to compare p6p to our result Eq. p2p above, we need to convert a general multipolar expansion of the Gauss 
curvature, say, 

RlK^l + Y,ClPl{^i) (37) 

l>2 

into a corresponding "height" expansion, say (denoting e; = {SR/Rq)i) 

^-E^'^^(M)- (3B) 

This is naturally done by defining e; so that the Gauss curvature on the 2-surface X(0, cj)) = {Rq + S R{^))'N {6 , (p) (with 
= sin0cos(/), A^^ = sin sine/), = cos6') embedded in an auxiliary 3-dimensional Euclidean space, coincides 
with Eq. (j37p . A straightforward calculation shows that this coincidence is obtained (in first order in the perturbation 
away from a round sphere) if 

ci = + 1) - 2]e, = {l- + 2)ei. (39) 



Note in passing that any eventual / ~ 1 contribution ei in psp is canceled in the corresponding curvature expansion. 

This agrees with the general result (j26bp above, and allows us to consider only multipolar orders / > 2. 

Finally, by combining Eqs. ([55)1 . . ([35]) and ([55)) . we find that the /-th shape Love number of a BH is given by 

In particular, we find for the first values of I: 

/i2 = -, h-i = — , /l4 = — . (41) 
4' 20' 84 ^ ' 

We see that hi diminishes rapidly as / increases. 

Note that all these Love numbers are smaller than 1. This contrasts with the Love numbers of a perfect- 
fluid star in Newtonian gravity which are given in terms of their "second" Love numbers ki by /i^cwton — ^ _|_ 2kf™^°'^ 
so that /iN°-«'ton > 1 (because fcf'™'™ > 0). In the recent work Q], it was, however, found that the shape Love 
numbers h^ of a general relativistic fluid star (say a neutron star) have a strong dependence on the self-gravity 
of the star, and that hf^ continuously decrease, as the "compactness parameter" c = GM/R increases, from the 
Newtonian value /jNcwton ^ ^ down to a value close to the above black-hole values, hf^, as the compactness reaches 
its maximum possible value. It was also found that the function hf^{c) becomes equal to hf^ as the compactness 
parameter c is formally continued to the black-hole value c^^ — 1/2. 

Using the Stirling approximation, one finds that, for large values of 



(42) 
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so that hi essentially decreases as 4 This rather fast decrease of hi, &sl increases, has useful practical consequences. 
Indeed, in the physically most relevant case where one is interested in the tides raised on a certain BH of mass M by 
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another BH (of mass m), located^ at a distance D, we see that the tidal expansion, Eq. of f/gxt (which, from a 
Newtonian point of view, would be expected to converge on the horizon as a geometric series ^^{Rq/ DY) raises on 
the horizon a tide whose height converges like 




which is essentially a geometric series ~ The appearance of ^ instead of implies a rather fast conver- 

gence, even when D is not much larger than Rq. Physically, this means that taking into account only the quadrupolar 
tide oc /i2 will probably suffice to give a good estimate of the full tide, even when the companion BH is rather near. 
In a later section, we shall, however, sec that this conclusion no longer holds when m <^ M and when the mass m is 
allowed to come very near the horizon of M. By contrast, when m ~ M, the finite size of the horizon of m docs not 
allow the distance D to become smaller than some minimum value I?min oc M + m, so that one might hope, in view 
of the remarks above, that the (linear) quadrupolar tide alone might give a good estimate of the tidal deformations of 
both holes down to the point where they formally touch. In the case of the quadrupolar tide, one knows from other 
studies m, [ill that a good measure of the / = 2 tidal field T2 is given by some component of the Riemann (or Weyl) 
tensor, which, in the exterior of a Schwarzschild BH of mass m is oc m/R^, where R is the usual (areal) Schwarzschild 
coordinate (rather, say, than the Weyl radial coordinate distance D that we have been using up to now.). This leads 
one to expect that, when m and M are comparable, m will raise on M a tide of height 

fSR\ u JJ^^p 1 m ( Ro{M) V 
\RoJ M ~ ^^AM\ R ) ' 

Upon "contact", i.e. when R = Rq{M) + i?o(m) = 2{M + to), we get {SR/Ro)m ~ j (Z+Cyi ■ This result would 
predict a maximum deformation SR/Rq 1/27, reached when M = 2m. We note that this deformation is surprisingly 
small. One should, however, remember that its "derivation" relied on several approximations that cease to be valid 
in the limit we consider. 

We have tried to further study the tidal deformations of comparable-mass BH's by considering the multi-black-hole 
Weyl solutions recalled in Eqs. ([8]) and (fTTj) above. The problem, however, is that these solutions do not allow one 
to study, in a physically relevant way, a process where two nearly isolated BH's get very near each other. Indeed, 
the total configuration will always include some struts or strings, that are necessary for global equilibrium. One can 
arrange two BH's, among N, to be free of such strings. However, when these two "central" BH's get close to each 
other, one will need other BH's (supported by strings) to become also very close to the central BH's. This generates 
additional strong forces and accelerations (related to / = 1 tidal terms) that mess up the pure I = 2 tidal heights 
coming from the (free fall) gravitational interaction of the two central holes. As a consequence, we could not use Weyl 
solutions to check the above prediction that "touching" BH's are only slightly deformed. 

Instead, we could, however, use well separated, multi-black-hole Weyl solutions to check the values of the first few 
Love numbers. For instance, one can take an asymmetric 3-black-hole Weyl solution, made of a central BH of mass 
M ("located" at z = 0, i.e. corresponding to a "rod" in the interval z € [—M,AI]), and of two "satellite" BH's: one 
of mass TOi located at z ~ bi > 0, and one of mass TO2 located at z ~ —62 < 0. One then considers a limit where 
61 ~ &2 tend to infinity (e.g. keeping toi ~ TO2 finite). In this limit, one can expand the solution ([8]) and (jlip in 
inverse powers of bi and 62. One finds that, by using a suitable constant C in Eq. (|8b[) . and by tuning TO2 (given mi, 
61 and 62) so that 

Toi o Toi TO9 n m2 „ / 1 \ , , 

one can arrange to have the central BH of mass M to hold in equilibrium without supporting struts (i.e. with 
limp^o 7 = on both sides of the z axis, and therefore with a smooth horizon and jd^^ = fi^^ = 1). Then, we find 
that the Gauss curvature of the central horizon reads 



^ Note the that extra Hnear terms n; 2Ti+2 j^'^f'^^ + ■•• that we neglected in Eq. I l29t are, in order of magnitude, essentially equivalent to 
correcting the distance D entering T; oc D~^~^ in the leading term riiTiM^ by D — » Z)(l C2 A/^/D^ -I- C4M'*/D*) + Such corrections 
include, in particular, possible effects of coordinate transformations on D. 
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Using the definition above of the Love numbers, one easily checks that this result is equivalent to /12 = 1/4 and 
/i3 = 1/20, in agreement with our general formula above. 



III. ELECTROMAGNETIC CASE: INFLUENCE OF FARAWAY CHARGES ON THE SURFACE 

CHARGE DENSITY OF A BH 

As a contrast to the gravitational polarizability properties of BH's, let us now consider their electric polarizability 
properties. For simplicity, we consider a physical situation where an uncharged BH of mass M is immersed in a 
general electric field generated by a static axisymmetric configuration of faraway charges. 
The background metric of the Schwarzschild black hole is 

ds^ = -{^- ^) '^^^ + ^ \m ^^^ + i'^^^ + ^^'^^) • (46) 

In the linearized approximation, the electromagnetic field F^^ = d^A^ — d^iA^ generated by the faraway charges must 
satisfy (outside the location of the charges) 

F^'r^ = 0, (47) 



where the semicolon denotes the covariant derivative. Eq. (|47|) yields a second-order partial differential equation 
for the scalar potential Aq = ~V. [Here V denotes the usual electric potential such that the electric field E = — 
V F = + V Aq.] Assuming that all the charges that generate Aq are at a distance D ^ 2M, we can, as in the 
gravitational case, consider that, in the intermediate domain 2M <C i? <C -D (where the spacetime is approximately 
flat), the potential Aq admits a flat-space multipolar expansion of the general type 

00 

A™'(i?,^) = ^T,i?'P, (cos 0). (48) 
/=o 

Here the coefficients r; ^ +d^Ao ^ +d^^^E arc the electric analogs of the tidal coefficients Ti ^ +d^U ^ +d^^^g 
that entered the asymptotic tidal expansion of the external gravitational potential J7ext , Eq. (jj]) . The sign convention 
is chosen so that a uniform electric field E = Egz directed along the positive z axis corresponds to ri = +E. In 
addition, the normalization of the t;'s is such that, in the particular case where the electric field that we consider is 
generated by a pointlike test charge q located at a large distance D along the z axis (in a flat region), r/ is simply 
equal to 

^^ = -^^ (49) 

SO that the asymptotic scalar potential reads 

00 

^ ^ - /,^',n a 75fe^'^'(-^^)- (50) 

yr'' + D'^ — 2rD cosO 

Let us now take into account the effect of the background spacetime curvature, Eq. Because of the spherical 

symmetry of the background, and the assumed axisymmetry, one can decompose the exact Aq = — y in a series of 
Legendre polynomials, say 



R 

1=0 1=0 



Ao = Y,ai{R)Pi{cose) = Y.\I'^' ^fwiiR)Pi{cose), (51) 



where the functions wi{R) = ai{R)/^l — 2M/R can be shown to obey the differential equation 

1 



(1 ^R^)^-^- 2R^ 
^ ' dR^ dR 



i-i?2 



wi = (52) 



where R = R/M — 1. Following [2^, we remark that Eq. ([52]) has the form of a general Legendre equation (m = 1), 
and admits two independent solutions, which can be expressed in terms of the associated Legendre functions of the 
first and second kind, Pl{R) and Q]{R), respectively. As a result, the expansion ([^T]) contains two classes of radial 
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functions ai{R): one class of solutions, say gi, is regular on the horizon, and behaves like i?' for R oo, while the 
other, say is logarithmically singular (except when ^ = 0) on the horizon and behaves like for R — * oo. We 

are interested here in the unique radial solution ai{R) that is regular on the horizon and grows like i?' when R ^ oo, 
i.e. the gi's, (so that the corresponding term in Eq. ((5T|) can match the term oc T;i?' in Eq. (|48p ). Such a solution is 
given by go = 1 when 1 = 0, and 



(20! 



dR 



when 1^0, 



(53) 



where the normalization of (|53|) is chosen such that gi{R) — i?' for R ^ oo. With this normalization, the unique, 
horizon-regular solution of (|T7|) which asymptotically agrees with pS]) (and corresponds, when / = 0, to an uncharged 
BH) reads 



(54) 



1=0 



Let us now consider what is the "charge density" cr [20| induced on the horizon of the BH by the electromagnetic 
field, namely 



1 

47r 



Decomposing cr into multipoles, 



dR 



R=2M 



1 

47r 



dV 
'dR 



a{cos9) = y^o-;P;(c 



1=1 



R=2M 



(55) 



(56) 



and using the identity 



one easily finds that (with Rq = 2M) 



ax 



1(1 + 1) 



x = l 



(57) 



(58) 



where 



uji + iy. 



I > 1. 



(59) 



On the left-hand side of Eq. ([55)1 we have introduced the quantity A-kR^cfi which has the dimension of an electric 
potential, i.e. the same dimension as the Z-th multipolar component T/i?g of the asymptotic electric potential 
(formally evaluated for a radius R = Rq = 2M). Eq. ((58)) is the electric analog of Eq. ((33)) : its right-hand side 
contains the i? — > i?o continuation of a "tidal" potential, while its left-hand side contains the "effect" of its tidal 
potential on the horizon (here the induction of a charge density). We can therefore consider that the dimensionless 
proportionality coefficient hf^ entering the linear relation l|58p is the electric analog of the hi Love number (hence 
our notation). 

It is interesting to note that this "electric Love number" is rather similar to its gravitational analog. It is related to 
it by the simple factor 



The first few values are hf^^ = 1, h 



EM 
2 



hf'" = 2{l-l)hi. 

1/2, h^^^ = 1/5. For large values of I, hf^^ decays as 



(60) 



EM 



lV2^l 



(61) 



As in the gravitational case, we can therefore expect that, as long as the inducing charge g is at a distance D > Rq, 
its "electric tidal" effect will be dominated by the lowest multipole, i.e., in the present case, by the induction of a 
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dipolar charge distribution Airai = hf^'n, i.e., (using our result hf^' ^ 1) Anai = i<;asympt .^^j^gj-g ^asympt ^j^g 
asymptotic {Rq <^ R <SC^ D) electric field strength. 

In this respect, it is also interesting to compare the electric Love numbers of a BH to those of a conducting 
sphere in flat space. We recall that a BH is analogous to a conducting sphere in that, at equilibrium, it is an 
equipotential surface [l^l, and that, in a general nonequilibrium situation, it exhibits dissipative properties (Ohm's 
law. Joule's law) similar to those of a conducting shell with surface resistivity equal to in = 0, H, 0, Q ■ The 

general "tidal" l-th multipolar component of the electric potential around an uncharged, conducting sphere of radius 
Rq (in flat space) is easily found to be 

Ao{R,9)^To + J2ji f^'-^j^Kcos^). (62) 

Computing the corresponding multipole-expanded charge density cr, Eq. (|55[) . we then deduce from the definition 
the flat-space values of the electric Love number of a conducting sphere 

^sphere ^ 2/ +1. (63) 

Note that these flat-space values are larger than 1 and they grow with I (ft,^P^°'^° ~ 3, h^^'^'^'^ ~ 5). This contrasts 
with the corresponding BH values {hf^ = 1, ^ ~ 1/2) that decrease with /. 



IV. ELECTROMAGNETIC CASE: NEAR-HORIZON LIMIT 



Up to now, we have been considering (both in the gravitational case and in the electromagnetic one) a BH 
immersed in the (gravitational or electromagnetic) "tidal field" generated by a configuration of faraway sources. To 
complete our understanding of tidal effects, it is, however, interesting to consider also the case where a (test) mass m, 
or a (test) charge 5, gets very close to the horizon. As a prelude to studying the gravitational case, we shall consider 
in this section the case where an external test electric charge get s very close to the horizon of a Schwarzschild BH. 
This was the situation studied long ago by Hanni and Ruffini [20] . 

We then consider a test charge q located^ at z = D along the positive z axis {9 = 0). Outside the radial 
location of the charge, the multipolar components ai{R) of Ao{R,d), Eq. ([5T|) . satisfy Eq. ((5^ (after the transforma- 
tion ai{R) = wi{R)y^l — 2M/R). Above, it was enough for our purpose to consider only the asymptotically-growing, 
horizon- regular solution of (|52p . for which a;(i?) is given by Eq. (|53p . Now, we need to consider also the 
asymptotically-decreasing, horizon-singular solution of (|52[) . for which ai{R) is given by [28j 

^'^^^ - -niTm^'^ ^ ^ — dR — ' ^^^^ 

The normalization of /; (i?) has been chosen so that 

fi{R) ^ R~^^+^\ R^oo. (65) 

In terms of the two radial solutions gi{R} and fi{R), one can write [2^ the electric scalar potential generated by a 
charge q located at z = as 

Ao = -V= -qY,[fi{D)gi{R)Q{D ^ R) + fi{R)gi{D)e{R - D)] Pi{cos0), (66) 
where Q{x) denotes the usual step function. 

Inserting the result ((66| into the definition ([SS]) of the induced charge density (and using again the result ((57|) used 
above to compute the electric Love numbers) leads to the following result for the l-th multipolar component of the 
charge density induced on the surface of the BH 

^^RoTi - -^^^MD)R'o = -hf"qfi{D)Rl (67) 



^ Note that, in this subsection, D denotes a Schwarzschild-coordinatc radial distance. 
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This result can be interpreted in two ways. On the one hand, one can view it as an appHcation of the general 
definition of Love numbers, Eq. (|58p . but with the understanding that, when the charge q gets near the horizon, 
one must replace the usual "EM tidal moment" r/ = —q/D^^^ by ti{D) = —qfi{D). On the other hand, one can 
alternatively (in the spirit of Love's original definition) wish to continue to define t; as being simply the "Coulombian" 
values Ti = —q/D^^^^ in which case one can say that the Love number h^^^ in Eq. (j58[) must be dressed by a distance- 
dependent correcting factor, say tf^^{D) = D^+\fi{D) (with ti{D) ^ 1 as £> ^ oo), so that hf^^{D) = hf^Hf^^{D) 
plays the role of a distance-dependent effective Love number in 

4.i?oa, = -hrhr\D)^R^, = -hr^iD)^Rl (68) 

Using Eq. ([M)) to evaluate fi{D) and remembering that [1^ Qi{x) ~ iP;(a;)log — P;_i(a;), where Pi_i{x) is 
a polynomial of order I — 1, one finds that the correcting factor tf^\D) = D^^^fi{D) grows, as D decreases from 
infinity down to _D = Rq = 2M , from 1 to a horizon value equal to 

Inserting the result ([M]) into the definition of the "dressed" , distance-dependent Love number leads to the following 
horizon value for the dressed Love number 

hf^^2M) = hf^Hf^\2M) = 21 + 1. (70) 

Note that, while the bare "faraway" Love numbers hf^'^ decreased roughly as 4~' as I increased, the dressed, near- 
horizon effective Love numbers (|70p now increase linearly in I. This linear increase with I is the multipole-expanded 
version of the result found in Hanni and RufSni that, as the charge q (with, say, q > 0) gets near the horizon, the 
induced charge density becomes nearly uniformly distributed on the horizon, apart from a strongly negative charge 
distribution, localized in a small patch "below" the charge q. Indeed, the limiting result ([70]) corresponds, when 
inserted in Eq. (j68p . to a charge distribution on the horizon given by (we recall that, by construction, the monopolar 
component ctq vanishes because the BH is uncharged) 

oo oo 

a = J2mPi{cose) = J2(^J + l)Pi{cos9) = [25(1 - fi) - I] , (71) 



1=1 ^' 1=1 







where (5(1 — fi) denotes a Dirac-delta distribution localized at = 1~ (with the convention d/x5(l — /i) = 1, Va < 1). 
Here we have (formally) applied the general theorem on the Lcgendrc expansion of a function f{fJ.) on the interval 
11 G [—1,-1-1], namely 

oo 

f{^i)=Y,fiPi{f^). (72a) 

1=0 

2/ + 1 f^^ 

fi = ^ / d^Kf{l^)m) (72b) 



to the distribution /(/^) = 25(1 — ^). 

Note also that the result (fZOll coincides with the Love number of a conducting sphere in flat space, Eq. ([63]) . This 
can be interpreted as meaning that, in the near- horizon limit, the electromagnetic interaction between the charge q 
and the horizon becomes localized in a small patch below g, and that such a localized behavior applies also to the 
interaction between a charge in flat space which becomes very close to a conducting sphere. 

V. GRAVITATIONAL CASE: NEAR-HORIZON LIMIT 

After this incursion into the electromagnetic analogs of tidal effects, let us come back to the gravitational case. We 
wish now to consider the gravitational analog of the Hanni- RufRni study [20| . i.e. the tidal deformation of the horizon 
by a test mass, m <^ M , in the limit where m gets very close to the horizon. To do this, it is convenient to start from 
the general result derived above for the curvature of the horizon, namely 



1 



{2MuYK{^i) = ^-d'^ . (73) 
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Correlatively, the cocfRcicnts q of the multipolar expansion of K{fi), 



{2MufK{^l) = 1 + ^c/Pz (cos 0), (74) 

1=2 



(where we used our general results (|26ap and (j26bp above) are given by 

21 + 1 



ci 



d^lPi{n)dl (1 - ^2-)g2(C/,.t(A.)-«) _ (75) 



Here, C/cxt(/^) is the value on the horizon of the Weyl-Newton external potential C/cxt, which is a solution of the 
axisymmetric, flat-space Laplace equation. We have in mind here a general configuration where t/gxt is generated by 
a stationary axisymmetric configuration of masses, m^, which include, among other sources, a test mass m <C M very 
close to the horizon. 

Let us consider the limiting case where all the external masses are test masses: rrii/M — > 0. In this limiting situation, 
C/cxt ( and u = = U^^) formally tends to zero, and we can replace Eqs. ([75)1 and (|75|l by their linearized 

approximations: 

(2M„)2i^"" - 1 = -dl [(1 - fi^) {tJ^M - ")] , ^ > 2, (76a) 
dfiP,{f,)dl[{l-ti')UMt^)]. (76b) 

In Eq. (j76bp we have used the fact that, for I > 2, the //-independent term u = U^l' does not contribute to q. 
Finally, as c|'", Eq. (|76b)) . is a linear function of [/oxt , we can consider that each c|'" is given by a sum over the various 
external masses, say, cj'" — cj'"(m) + c\"^{mi). In the following, we shall focus on the contribution cj'"(m), where 
TO is a test mass which is very close to the horizon ^. The individual contribution c}™(to) is obtained by inserting in 
the right-hand side of Eq. (j76bp the p limit of —ipm, where tprn is given by Eq. (|10ap . As we are considering the 
test- mass limit m/M 0, we can simplify this logarithmic expression for = ~V'm by expanding the logarithm 
in powers of 2m/{R^ + R^), and then by approximating (i?^j + i?~j)/2 simply by p'^ + (z — b)^ . This yields the 
simple, Newton-like result 

~ ^ (77) 

Then, taking the horizon limit (i.e. p 0) and replacing z by Mp according to Eq. (|16p leads to the explicit result 
(we assume, for definiteness, that h > M > 0, i.e. that the mass to is near the north pole, p = +1): 



91 + 1 

c'rim) ^ / dpP,{p)d'^ 



to(1 — p'^) 
h-Mp 



(78) 



There are now two ways to discuss what happens when h decreases from h ^ M down to ~ Af . [We recall that h 
is the Weyl-coordinate distance, so that the horizon is located at 6 = M ] A first way consists in noticing that the 
integral ([78|) can be explicitly expressed in terms of the Legendre functions of the second kind. Indeed, using [2^ 

Qi{^) = \ r Piiv)^, (79) 



cf-{m)^{2l + l)^{b'~l)dlQi{b), (80) 



we can reexpress Eq. ((78| as 



where b = b/M. 

Similarly to our discussion of the electromagnetic case, we can then re-interpret the result (|80[) by writing the l-th 



^ Though the present linear approximation allows us to focus, independently from the other masses rrii, on the contribution cj'"(m), 
one should keep in mind that the various masses are not totally independent from each other as they must respect the "equilibrium 
condition" Jsj, which is needed for ensuring the regularity of the metric on the portions of the z axis that touch the "central" BH of 
mass M. 
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il-l)il + 2)hMb)[^^] =i2l + l)ib'-l)diQiib). (82) 



multipolar component of the adimensionalizcd horizon curvature "raised" by an external test mass m at Weyl distance 
b in the form 

[(2M,)^X""], c cf"(m)PKM) = {l- m + 2)^'*^^) (^:^^yTT^^'(/^)' (81) 

where hi is the bare, faraway gravitational Love number, Eq. (|34[) . and where ti{b) is a distance-dependent correcting 
factor (normalized so that ti{b) ^ 1 as 6 ^ +00): here we conventionally replaced the factor in the electro- 

magnetic definition by l/{b + A/)'+^ to ensure that these "Coulombian" factors agree both when b +00 and 
on the horizon. 

Comparing Eq. ^ to Eq. we get 

{l-l){l + 2)hiti{b) f- 

If we now formally^ \ct € = b — M tend to zero, i.e. 6 = b/M = 1 + e (with e = e/M) tend to 1, one finds that 

ip - \)dlQi{b) ^ i [1 - Z(/ + 1)6 + Oil")] , (83) 

the crucial point being that the coefficient of the leading term = (6 — 1)~^ in the right-hand side is equal to 1, 
independently from the value oil. In other words, the near- horizon limit of the multipolar coefficients of the curvature 

is (e = - M) 

cf"(m)|b^M = (2/ + 1)-, (84) 
e 

and this result can be interpreted (similarly to the electromagnetic case (|70p ) as coming from the horizon limit of a 
distance-dependent "dressed" Love number hi{b), 

21 + 1 M 

HM^e)^hiriiM^e)^^^-^^-. (85) 

In other words, the correcting factor nib) grows, as 6 ^ M , in a strongly ^-dependent manner (roughly ^ 4+'), and 
this growth compensates the usual behavior [hi ~ 4~') of the faraway Love number. The second way to obtain the 
simple final result is to explicate the second derivative entering Eq. ([75]) so as to obtain: 



cl-(m) = (2/ + l)^ f'dnPiitJi) 



52-1 



(86) 

(21 + 1)- dpiPiifi)- 



e y_i *^-'(l + 6-/.)3 

Then one notices that, in the near-horizon limit e = e/AI 0, the function 5i{l — ^) ~ 26^(1 + e — = 
9p[e2(l-|-e— /i)~2] gets localized near /i = 1^ and has an integral diiSi{l — fj.) = [i^{l+i—iJ.)~^]l^ = l — e^ll+e—a)'^ 
which tends to 1. In other words. 5i{l — /i) 5(1 — fi) as e — * 0. Using Pi{l) = 1, we then obtain another proof of 
the limiting result (|84)) . 

Inserting the result ((84|) in the multipolar expansion ([74|) . we then see that the contribution to K coming from a mass 
m tends, as m approaches the horizon (e = b — M <^ M) towards 

oc 

(2M„)2ii'i-(m) ~ - V(2/ + l)Pl{^i). (87) 

1=2 

Using, as we did in the near-horizon electromagnetic case above, the multipolar decomposition of the Dirac-delta 
distribution 25(1 — /^), we can sum the scries ((57)) to get 

(2Af„)2i^'-(m)^^^^^ [25(1 (88) 



* As our analysis here assumes f/ext <S 1, we should always keep m/e <§; 1. In other words, the coordinate distance to the horizon, 
e = b — M, should tend to zero less rapidly than m/M. 
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This result is the gravitational analog of the electromagnetic result ([71]) . It shows that, modulo a global smooth 
curvature linear in /x, the tidal horizon curvature generated by a near-horizon mass m is localized in a small patch 
"below" the mass m. Similarly to the electromagnetic case, this result could also be formulated in saying that, in the 
near-horizon limit, the tidal curvature, being localized in a small patch on the horizon (as witnessed by the integrand 
oc (6^ — l)/(6— jj,)^ in Eq. ([5^ ). could also be computed by replacing the real curved horizon by a flat-space horizon 
or "Rindler horizon" (as done in [1, HH). 

We will, however, get a more interesting result (which has no electromagnetic analog) by considering the "shape 
distortion" of the horizon associated to the tidal curvature (j88| . Using our general result ([39|) . the multipolar 
expansion of the corresponding near-horizon shape distortion is 

6R{^i) m ^ 21 + 1 



R, e ^ /(/ 

To sum this new series, we apply the differential operator A -I- 2, where A denotes the Laplacian on the unit sphere 
(with AP/(/i) = —1(1 + l)Pi{fi)) to Eq. (15^ and find that SR{fi) satisfies the differential equation 

-(A + 2)^ = l + 3Ai-25(/i-l), (90) 

TO Rq 

where A, acting on an axisymmetric function, can be replaced by 9^(1 — /i^)9p. This yields a second-order ordinary 
differential equation for 6R{ii). The presence of a delta- function singularity on the right-hand side indicates that 
dR{fi) must have a ~ ln(l — fj,) singularity at /i 1 (and no singularity at /i — > —1, contrarily to the Qo{fJ.) Legendre 
function, which satisfies (A -I- 2)Qo = 0, modulo two delta-singularities at /i = 1 and fi = —1). One then easily finds a 
particular solution of ([90|l of the form (A + 2) (— ^ — /iln(l — /x)) = 1 -|- 3^ — 2(5(/i — 1). The general solution of (^0]) 
is then obtained by adding to this particular solution a general regular solution of (A + 2)/(/i) = 0, i.e. a multiple of 
Pi{^i) = fJ-. One can determine this multiple by requiring (as needed from the absence oi I = and I = 1 contributions 
|)) that j diJLPi{^)5R{ii) = 0. Finally, this uniquely determines the sum ([89]) to be 



SR(p) m 
Ro ^ ~ 



i-Ailn(l-Ai)+Ai(ln2-^) 



(91) 



The corresponding shape is illustrated in Fig. 1 (for a prefactor m/e = 1/30). This Figure shows the "height" 
of the tide raised on the horizon by a mass m, in the limit where the location b of the mass tends to the horizon 
{b — > M, keeping m/e = m/{b — M) finite and small). When e is kept non zero (even beyond the prefactor m/e), 
one finds that the logarithmic "spike" at = cos6' = 1 (north pole) is rounded off on the characteristic angular scale 
6c ~ \/ {b — AI)/M, corresponding to the multipolar orders Ic ^ ^ ^/Alje. [The characteristic multipolar order 
Ic can be read off the ncxt-to-leading term in the expansion on the right-hand-side of Eq. (|83p .] 



VI. CONCLUSION 



We compared and contrasted the gravitational polarizability properties of black holes (BH) with their electromagnetic 
polarizability properties. Our main results are: 

1. The definition and computation of the infinite sequence of "shape" (or "height") Love numbers hi of a BH, 
where I = 2,3,4... is the multipolar order. The result is given by Eq. (|40|) : hi essentially measures the ratio 
between the "height" of the l-th tidal bulge raised on the horizon of a BH and the corresponding external l-th 
tidal potential C/oxt(Po) (analytically continued down to the horizon Rq ~ 2GM). Contrary to the Newtonian 
hi Love numbers of a perfect-fluid star which are larger than 1 0, we found that the Love numbers of a BH 
are smaller than 1, and tend rapidly (and exponentially) toward zero as I increases, e.g. /12 = 1/4, /13 = 1/20, 
/14 = 1/84,..., hi ~ \/27rI/2^'+^. In a related recent work fj^, it was found that the height Love numbers of a 
neutron star approach those of a BH as the "compactness" c = GM/R of the star formally tends toward the 
BH value c^" = 1/2. 

2. The corresponding definition of a sequence of electromagnetic Love numbers hf^^ {I = 1, 2, 3, ...). These essen- 
tially measure the ratio between the l-th multipolar charge induced on the horizon of a BH and the corresponding 
external l-th multipolar electric potential (analytically continued down to the horizon). Here again we found 
that the electric Love numbers of a BH (given in Eq. ([59]) ) are much smaller than those of a conducting sphere 
in fiat space ((SS]) . and tend rapidly (and exponentially) toward zero as I increases: hf^^ = 1, ft.f'^^ = 1/2, 
hf^^ = 1/5, hf^'^ ~ /-\/27r//2^'. In addition, we found that they arc simply related to the gravitational 
height Love numbers: hf^' = 2{l — l)hi. 
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FIG. 1: The shape of a meridian section of the "tidal bulge" raised by a test-mass m located at the north pole of the horizon 
of a BH. 



3. The comparative study of the gravitational and electromagnetic polarizability properties as the "polarizing" 
mass or charge approaches the horizon. Both cases can be described by replacing the "bare", faraway Love 
numbers, by some "dressed", distance-dependent factor ti{b)hi or tf^^{D)hf^^. We found that the gravitational 
(respectively, electromagnetic) dressing factors ti{b) (respectively, tf^\D)) compensate, in the near-horizon 
limit, the exponential decrease of hi (resp. hf^^), and lead, in the gravitational case, to a specific "spiky" shape 
for the distorted horizon given by Eq. ([9T|) . and illustrated in Fig. 1. 
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